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An ingress and a complete transit of HD 80606 b 
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ABSTRACT 

We have used four telescopes at different longitudes to obtain near-continuous light 
curve coverage of the star HD 80606 as it was transited by its ~4-Mj Mp planet. The 
observations were performed during the predicted transit windows around the 25 th 
of October 2008 and the 14 th of February 2009. Our data set is unique in that it 
simultaneously constrains the duration of the transit and the planet's period. Our 
Markov-Chain Monte Carlo analysis of the light curves, combined with constraints 
from radial-velocity data, yields system parameters consistent with previously re- 
ported values. We find a planet-to-star radius ratio marginally smaller than previously 
reported, corresponding to a planet radius of R p = 0.921 ± 0.036i?j„ p . 
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1 INTRODUCTION 

An important class of the ~ 400 extra-solar planets known 
to date are the so-called "hot Jupiters", orbiting a fraction 
of an AU from their host star. Many of these have signifi- 
cantly larger radii than planets of similar mass in our Solar 
System. This has been a challenge for theoretical models of 
their structure. The intense radiation and tidal forces from 
the host star are expected to play a signific ant role (e.g. 
iBodenheimer et alJliuOll ; iBurrows et al.ll2007l V For planets 
in highly eccentric orbits, both of these factors vary greatly 
over the course of the orbit. Thus, such planets provide in- 
teresting tests for models o f the structure an d dynamics of 
planetary atmospheres (e.g. llrwin e t al. 2008i, and references 
therein). 

A planet in a 111-day orbit around the star HD 80606 
was fi rst detected in radial velocity observations (|Naef et al.l 
l200ll h The minimum mass (M p sini) of the planet is 3.9 
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times the mass of Jupiter. HD 80606 b has the most ec- 
centric orbit of all the extra-solar planets known to date 
(e = 0.93). Infrared (8/im) observations clearly show the 
rapid he ating of the planet's a tmosphere during periastron 
passage (|Laughlin et al.l 120091 . herein referred to as L09). 
Laughlin et al. also reported detection of a secondary eclipse, 
implying the orbital inclination is close to 90 degrees, and 
motivating efforts to observe a transit of the planet in front 
of the star. 

In early 2009, several grou ps observed a tran- 
sit egress in photometry (|Moutou et al] |200£ ; 
iGarcia-Melendo fc McCulloughl 120091: iFossev et all I200S ) 
and spectroscopy ( Moutou et all 12009). Analyses of these 
data, together with old and new radial-velocity mea- 
surements provided the first constraints on the planet's 
radius and actua l mass, as well as its orbita l parameters 
(|Pont et all 120091 . herein referred to as P09; iGillonl 120091 . 
G09). These analyses are limited by the fact that the dura- 
tion of the transit is not constrained, whi ch in turn in c reases 
uncertainties in the system parameters. IWinn et al. I (|2009l . 
W09) report multiple observations of the transit ingress in 
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June 2009, and combine their data with the earlier ingress 
observations and Keck radial-velocity data to constrain the 
transit duration. 

In this paper we present previously unpublished multi- 
site photometric observations of a complete transit of 
HD 80606 b on 14 February 2009, and an ingress on 25 Oc- 
tober 2008. We describe the observations and the resulting 
data sets in section [2] In section [3] we describe the methods 
we used to analyse our data, and in section [4] we report and 
discuss the results. 



2 OBSERVATIONS 

We observed HD 80606 with four telescopes, during pre- 
dicted transit windows around 25 October 2008 (herein re- 
ferred to as Oct08), and 14 February 2009 (Feb09)Q As no 
transit or secondary eclipse had been reported at the time, 
these windows were ~ 24 hours long. 

Where available, we observed through Bessell B and R 
filters, alternating between them every few minutes. The one 
exception to this is noted below. In all cases we deliberately 
defocused the telescope (often to the point where the point- 
spread function became an annulus, but taking care that the 
target and its binary companion do not blend) . This allowed 
longer exposure times without saturating the detector. More 
importantly, it meant that each star was measured by a large 
number of pixels, reducing the effects of flat-fielding errors 
and pointing drift on the photometry. 

Besides those mentioned below, we also attempted ob- 
servations at Bisei Astronomical Observatory in Japan, Kry- 
onerion Astronomical Station in Greece, and an amateur 
observatory in Athens, Greece, but did not yield any useful 
data due to adverse weather. 



2.1 Faulkes Telescope North, Hawaii 

The 2 m Faulkes Telescope North (FTN) — part of 
the Las Cumbres Observatory Global Telescope Network 
(LCOGT.net) — is located on Haleakala in Hawaii. We 
used FTN and a Merope camera with a 2kx2k e2v CCD 
(binned 2x2). Observations were conducted by the tele- 
scope's Robotic Control System, alternating between sets of 
37-sec exposures in B and 20-sec exposures in R. We only 
obtained useful data for the Feb09 transit. Images were pro- 
cessed by our automated pipeline to remove the effects of 
bias, dark current and CCD sensitivity variations. 



2.2 IUCAA Girawali Observatory, India 

Observations were carried out from IUCAA Girawali Obser- 
vatory (IGO; lDas et~ai]|l999h and we obtained useful data 
during the Feb09 transit window. 

We used the IUCAA Faint Object Spectrograph Cam- 
era (|Guptall2002j . IFOSC) mounted on the direct Cassegrain 
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focus of the IGO 2 m telescope. Exposure times were 15- 
30 sec in B and 4-10 sec in R. Images were corrected us- 
ing overscan strips and master bias frames followed by flat- 
fielding, using the IFOSC photometry pipeline package de- 
veloped under IRAF. 



2.3 Oskar-Luhning Telescope, Germany 

We observed HD 80606 on three nights in Feb09 with the 
1.2 m Oskar-Luhning Telescope (OLT) at Hamburg Obser- 
vatory. We changed filters after every 5 exposures, with ex- 
posure times of 55 s in B and 35 s in R. Images were bias- 
/dark-subtracted and flat-fielded using the Figaro subrou- 
tine in the Starhnk software package. 



2.4 Liverpool Telescope, Canary Islands 

Liverpool Telescope (LT; Isteele et alj|2004h data were ob- 
tained on 2008 October 2 5 and 26 using t he RISE frame 
transfer wide field camera (|Steele et alj [2008) in 2 x 2 binned 
mode, giving a pixel scale of 1.07 arcsec/pixel. The fil- 
ter employed consisted of 2 mm Schott KG5 bonded to 
3 mm Schott OG515, giving half-maximum transmission 
from 515 nm to 695 nm (i.e. roughly equivalent to a com- 
bined V + R filter). An exposure time of 1 sec was used, 
with effectively no dead time between exposures. All data 
were dark-subtracted and flat-fielded using the ULTRA- 
CAM pipeline, which is optimised for time series photome- 
try. A master twilight flat was constructed from a median 
of several hundred frames. 



2.5 Photometry 

Starting with the calibrated images, we performed aperture 
photometry on all data sets in the same manner, with the 
exception of the LT data. 

For each instrument, a single master image was se- 
lected, and used as a reference for all images obtained with 
that instrument. The coordinates of the target and refer- 
ence stars in the master image were obtained by convolving 
it with a copy of the PSF from the same image, then running 
DAOFIND in IRAF on the result. 

All other images were registered to the master using 
direct correlation with the XREGISTER task. The master 
list of coordinates was then transformed into the reference 
frame of each image. DAOPHOT was used to obtain pho- 
tometry in a range of apertures (slightly different for each 
instrument). The photometric scatter was smallest in the 
largest aperture used for all data sets (radius 11" for FTN, 
8" for IGO, and 10" for OLT). 

For the LT data, we also performed aperture photom- 
etry, trying out different aperture radii, choosing the one 
that minimised the noise (8.6 arcsec). Apertures were cen- 
tred in each frame by cross-correlation of the image with a 
Gaussian. 

We measured the flux of HD 80606 relative to its binary 
companion HD 80607, as it was the only reference star avail- 
able in all data sets. Additional stars in the fields of view 
of some of the instruments were significantly fainter and of 
different colours, and were only used as check stars. 
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Figure 1. The complete data set analysed in this paper: A transit ingress of HD 80706 b observed with the Liverpool telescope (LT, 
top panel), and the following transit observed by the Faulkes Telescope North (FTN), the IUCAA Girawali Observatory's 2m telescope 
(IGO), and the Oskar-Liihning Telescope (OLT). The middle and bottom panels show the same transit in the Bessell R and B bands, 
respectively. The horizontal bar across the top shows the section of data used in the fit. The error bars indicate the weights used in the 
fit, i.e. they include the statistical uncertainty as well as the red noise estimate (pi and oy iS in equation [TJ . 



We prepared each light curve for fitting by taking the 
following steps. 

(i) We removed points more than 5a from the median of 
the time series. 

(ii) We normalised the flux values to a median of 1 for 
points well outside the transit. 

(iii) We binned the light curve into 10-minute bins in or- 
der to reduce the number of points to fit. We estimated 
uncertainties as the standard deviation of the values within 
each bin, divided by the square root of the number of points 
in the binQ 



3 ANALYSIS 

The complete set of useful data we obtained is shown in 
Figure [1] It includes a transit ingress observed by the LT in 



2 For bins containing only a single original point, we estimated 
the uncertainty based on the nearest two bins with at least two 
points. 



October 2008, and — between the FTN, IGO and OLT data 
sets — nearly continuous coverage of the transit in February 
2009 in two filters. 



3.1 Light-curve model 



We modelled the HD 80606 system as a planet of radius R p , 
orbiting a star of radius R t , with orbital period P, semi- 
major axis a, eccentricity e, and inclination i (relative to 
our line of sight). The argument of periastron is u, and the 
planet passes there at time T p . 

To obtain light curves for a given set of system pa- 
rameters, we first solved Kepler's equation using Newton- 
Raphson iteration, then calculated the planet's position in 
its orbit, and its distance (in the plane of the sky) from 
the centre of the star. Taking thi s as in put, the analytical 
transit models of iMandel fc Agoll (|2002T ) provided the rela- 
tive fluxes expected for our physical model. We assumed a 
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quadratic limb-da rkening law, w ith coefficients interpolated 
from the tables of lClaretl i fcOOcM 5 ! 

3.2 Fitting 

We determined the parameters of the model using a Markov- 
Chain Mont e Carl o (MCMC) approach similar to that of 
I Torres et al.l (|200&f ) and references therein. Briefly, from a 
random starting point, small jumps in parameter space are 
generated, and each new point is evaluated using a merit 
function (x 2 )- If the jump results in a lower x 2 > it is ex- 
ecuted (i.e. the new point becomes the next "link" in the 
chain). Otherwise, it is only executed with a probability 
exp(— Ax 2 /2). If a jump is not executed, the previous point 
is repeated in the chain. After a large number of jumps, the 
result is a distribution of points in parameter space that ap- 
proximates the joint posterior probability density of all the 
parameters given our data. 

We used a merit function of the form: 

2 _ ijlsfj - mf ^ (Vj - Cjf 

The first sum is the usual chi-squared statistic and is 
related to the probability of the data given the model. For 
each of iV data points, /, is the measured relative flux, rrn 
is the model at that point, Oi is the statistical uncertainty 
in the measurement, and oy, s is an estimate of the corre- 
lated noise in each data set (see below). We also included a 
normalisation factor n s (constant within each data set s) to 
allow the model to match the out-of-transit relative flux. 

The additional terms in the merit function represent 
constraints from previously published measurements, as de- 
tailed in Table [1] In each case a value calculated from our 
model (vj) is being compared to a measured value (cj) with 
corresponding uncertainty (o c ,j)- We included the time and 
duration of the secondary eclipse from L09, as well as the pe- 
riod, eccentricity and argument of periastron they derived 
from Keck radial-velocity data. W09 note that the uncer- 
tainty in uj quoted by L09 is likely to be underestimated. 
Indeed, both W09 and P09 analyse larger data sets and find 
uncertainties in uj of ~ 0.2°. Therefore we adopt the value 
of uj from L09, but only give it a weight corresponding to 
Oc.u; = 0.2. We did not include constraints from any other 
photometric observations in order to obtain results as inde- 
pendent as possible from other analyses. 

During the MCMC simulation we varied the following 
parameters: a/R t , P, T p , e' (see below), uj, b = a/R*cosi 
(proportional to the transit impact parameter), R p /R*, and 
a normalisation factor for each instrument/filter combina- 
tion (n s ). A total of 14 parameters were varied. In order 
to avoid strong correlations which restrict the movement of 
the MCMC chain through parameter space, we used com- 
binations of some of the physical parameters that are more 
directly constrained by the data (e.g. the radius ratio in- 
stead of the absolute radii). When eccentricity was included 

3 Based on the stellar parameters reported bv lNaef et al. 
we obtained the limb-darkening coefficients «i = 0.5043, 112 = 
0.2535 in R band, and m = 0.8827, u 2 = -0.0166 in B band. For 
the LT data set (V + R), we used the mean of the coefficients for 
V and R: m = 0.6117, u 2 = 0.1853. 



Table 1. Measurements by lLaughlin et al. I [120091) included as 
constraints in the MCMC merit function. The three columns cor- 
respond to Vj, Cj and cr c j, respectively, in equation [T] The last 
two rows refer to the epoch and duration (2 nd to 3 rd contact) of 
the secondary eclipse. 



Parameter 


Value 


Uncertainty 


Orbital period, P (days) 


111.4277 


0.0032 


Orbital eccentricity, e 


0.9327 


0.0023 


Argument of periastron, uj (deg) 


300.5 


0.2° 


Epoch of mid-eclipse (HJD) 


2454424.736 


0.003 


Total-eclipse duration (days) 


0.070 


0.009 



a This is our conservative estimate of the uncert a inty (see 
section 13.21 1. The value reported by lLaughlin et al. (2009) 
uj= 300.4977 ±0.0045°. 



as one of the varied parameters, correlations caused differ- 
ent MCMC chains to converge to different locations in (P, 
e, uj) space. We found that most of the variation in e over 
all chains could be described as a linear function of Pand 
uj. Thus, we replaced e with e' = e — koP — kiuj, where 
fco = 0.035 and fci = —0.0035 were determined from the 
previous fit. This resulted in the MCMC simulation moving 
more freely through parameter space and a more continuous 
distribution of points when the chains were combined. 

We generated jumps by adding a Gaussian-distributed 
random value to one of the varied parameters (randomly 
chosen) at each iteration. We tuned the jump sizes for each 
parameter so that approximately 25-50% of jumps were ex- 
ecuted. 

3.3 Correlated noise 

Each of our data sets contains some noise that is correlated 
in time due to gradual variations in atmospheric or instru- 
mental factors. The internal uncertainties we estimated for 
each data point (as described in section |2~5")) do not measure 
this red nois e, yet it can sign ificantly affect the results of 
our analysis (|Pont et al.| [2006) . To account for the effect of 
red noise, we reduced the weight of each point in the merit 
function (equation [TJ with an additional noise term, oy jS , 
which we estimated as follows. We first performed an initial 
fit setting oy jS to an optimistic 0.0005 for all data sets. For 
each data set s (of N s points), we then found the value of 
ay, s for which the reduced chi-squared 

i — l 

was equal to unity0 We then repeated the fit, and checked 
that the values of oy jS were self-consistent. The final 
values used were LT: 0.001, FTN(P): 0.0016, IGO(P): 
0.0015, OLT(P): 0.0012, FTN(B): 0.0019, IGO(P): 0.0017, 
OLT(B): 0.001. The average internal uncertainty of data 
used in the fit was only 0.0006, so — as is generally the case 
in such analyses — red noise is the main factor limiting the 
accuracy of our results. 

4 The effective number of parameters for this chi-squared is 8, 
consisting of the 7 physical parameters and one normalisation 
factor n s for the given data set. 
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Figure 2. Close-up of the data used in the fit, along with the best-fit model obtained from MCMC simulations. Symbols and error bars 
are the same as in Figure [l] only the horizontal scale is different. 



Table 2. Physical parameters of the HD 80606 b system obtained from the MCMC fit. 



Parameter 


Median 


Uncertainty 


Scaled semi-major axis, a/R* 


101.2 


-5.5, +6.9 


Orbital period, P (days) 


111.4273 


±0.0031 


Orbital eccentricity, e 


0.93369 


±0.00068 


Epoch of pcriastron, T p (HJD) 


2454424.8575 


±0.0040 


Argument of periastron, ui (deg) 


300.53 


±0.19 


Orbital inclination, i (deg) 


89.341 


-0.063, +0.073 


Planet/star radius ratio, R p /R t 


0.0967 


-0.0035, +0.0032 


Planet radius, R p a (Rj U p) 


0.921 


±0.036 


Epoch of mid- transit (HJD) 

Complete transit duration ((4 — t\ , hours) 

Ingress/egress duration (i2 —t\, hours) 


2454876.3173 
12.14 
2.33 


±0.0036 
±0.36 
±0.37 



a Using R, = 0.978 ± O.O15i? l|Pont et al.ll2009l V 



4 RESULTS AND DISCUSSION rameter we report the median of this distribution as the 

For our final fit, we created 25 independent chains of 42,000 
steps each, and discarded the first 2,000 steps of each chain 
to minimise the effect of the (randomly chosen) starting 
points on our results. Combining the remaining steps of all 
the chains, the resulting distribution of parameter values 
approximates their joint probability density. For each pa- 



6 Hidas et al. 



"best-fit" value0 and the boundaries of the central 68% as 
the l-cr uncertainties. 

The results thus obtained are shown in Table [2] and 
the fit itself is compared to the data in Figure [2] The re- 
duced chi-squared for this fit is 0.63, which suggests we have 
slightly overestimated the uncertainties in our data. Thus 
the uncertainties in our results are somewhat conservative. 
We also note that there are unavoidable correlations between 
some parameters. Besides the P-e-ui correlation mentioned 
in section T3.2I the strongest of these are between a/R t , b 
and R p /R*. There are also significant correlations between 
T p and uj, and also between R p /R* and the normalisation 
factor for the IGO data. The uncertainties in results include 
the effects of all these correlations. 

Since we used the values of the period (P), eccentricity 
(e) and argument of periastron (uj) from L09 as direct con- 
straints in the fit, we obtained results that agree with those 
values. We did find a larger uncertainty for w, but a more 
precise value for the eccentricity. 

Our results also agree well with the parameters reported 
by P09, G09 and W09. The only significant exception to this 
is the period, where our result is 3a smaller than the other 
analyses. We note, however, that we have essentially forced 
our simulations to agree with the period reported by L09. If 
we repeat our fit without the direct constraints on P, e and 
uj (see section [372|) . the best-fit period is almost unchanged, 
but the uncertainty in its value is four times greater. This 
means our photometric data are in fact consistent with the 
period reported by P09, G09 and W09. 

We also find a planet-to-star radius ratio that is 
marginally smaller (by 1.8cr) than the values reported by 
W09 and G09. A possible cause for this is the nature of 
our data from the IGO, which contains the majority of 
the in-transit points, but no out-of-transit points on the 
same night. This means the normalisation of those data is 
not strongly constrained, and leads to the correlation noted 
above. In addition, the IGO data appear to show a signifi- 
cant systematic trend in the middle of the transit, which we 
were unable to remove by de-correlation against any physical 
parametrers of the observations (such as airmass or position 
on the CCD). 



5 CONCLUSION 

We have obtained multi-site observations of a transit ingress 
and a complete transit of HD 80606 b across its host star. 
We analysed these data independently of any other photo- 
metric data, and found system parameters consistent with 
previously reported values. 

These observations were made using four telescopes at 
different sites. This allowed us to obtain near-continuous 
coverage of this 12-hour event. However, the differences be- 
tween the instruments, telescopes and time-allocation proce- 
dures were, to an extent, limitations on our ability to obtain 
a uniform data set. In the near future, the completion of 
LCOGT's network of lm robotic telescopes (jBrown et al. I 

5 As lPont ct al. (2009) point out, the set of parameters thus ob- 
tained does not actually correspond to the MCMC step with the 
lowest x 2 value. This is because each parameter's median value 
comes from a different step in the chain. 



|2010| ) will greatly facilitate observations of this kind, pro- 
viding near-identical instrumentation at a number of sites, 
under the control of a flexible, central scheduling system. 
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